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PROCEDE D' ESTIMATION D'UN RAYONNEMENT DIFFUSE, NOTAMMENT AFIN DE CORRIGER DES 
MESURES EN RADIOGRAPHIE. 

La radiographie d'un objet est amelioree en estimant 
le rayonnement diffuse qu'il transmet aux detecteurs. Pour 
cela, on se sert du rayonnement diffuse effectivernent me- 
sur§ & travers un simulacre de I'objet, presentant des pro- 
priety d'attenuation analogues, et qu'on modifie par des 
coefficients de ponderation obtenus par une transformation 
des valeurs du rayonnement total repu & travers I'objet (3) et 
le simulacre (8) selectibnne. On parvient ainsi & ameliorer 
les radiographies sans faire subir une double irradiation & 
I'objet pour mesurer s£par6ment le rayonnement diffuse. 
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PROCEDE D' ESTIMATION D'UN RAYONNEMENT DIFFUSE , 
NOTAMMENT AFIN DE CORRIGER DES ME SURE S EN RADIOGRAPH IE 

DESCRIPTION 

Le sujet de cette invention est un procede 
d' estimation d'un rayonnement diffuse, dont 
1' application principalement envisagee est la 
correction de radiographies . 

L ' utilisation d'un rayonnement conique, 
tr6s fr^quente en radiographie, presente 1 ' inconvenient 
de produire un rayonnement diffuse important a travers 
1'objet examine- En d'autres termes, chacun des 
detecteurs situes derriere 1'objet regoit non seulement 
un rayonnement primaire, provenant directement de la 
source par un trajet rectiligne et ayant traverse une 
region bien definie de 1'objet, mais un rayonnement 
diffuse de provenance indeterminee qui affecte la 
mesure et qu'il serait done souhaitable de corriger. 

Plusieurs proc^des sont deja pratiques. 
C'est ainsi que .le rayonnement primaire peut §tre 
mesure seul si une collimation stricte des detecteurs 
et de la source est faite afin d' intercepter le 
rayonnement diffuse, mais ce procede necessite en 
pratique un balayage du faisceau qui est lent a 
accomplir, et pendant lequel on doit s ' accommoder de 
mouvements du patient si on examine des etres vivants . 

On a aussi eu l'idee contraire de ne 
mesurer que le rayonnement diffuse. On dispose pour 
cela un reseau discontinu d' absorbeurs , comme des 
billes de plomb, entre 1'objet et les detecteurs, pour 
arreter localement le rayonnement primaire, de sorte 
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que les detecteurs situes derriere ces absorbeurs ne 
mesurent que le rayonnement diffuse. Ce procede appele 
« beam stop » donne done des tables ou nappes 
bidimensionnelles de valeur de rayonnement diffuse, 
qu'on complete par interpolation entre les detecteurs 
places derriere les absorbeurs. Le rayonnement diffuse 
ainsi estime est soustrait du rayonnement total mesure 
separement. Ce procede est precis mais a 1 ' inconvenient 
qu'il impose deux irradiations de l'objet et done un 
doublement de la dose de rayons qu'il "recoit. Un 
dernier exemple de methode de correction du rayonnement 
diffuse par des moyens materiels comporte l'emploi de 
grilles anti-dif fusantes, mais leur efficacite n'est 
que partielle ; elle est insuffisante pour un faisceau 
conique, ou le rayonnement diffuse peut etre plusieurs 
fois superieur au rayonnement primaire. 

Enfin, il existe un certain nombre de 
methodes nurneriques pour estimer le rayonnement 
diffuse, a partir de convolutions ou de de convolutions 
des mesures par exemple ,- on pourrait aussi citer le 
brevet frangais 2 759 800 pour un procede numerique 
different, analytique. Elles sont en general d'emploi 
delicat car elles dependent de parametres choisis par 
l'utilisateur (noyaux de convolution par exemple) qui 
ne donnent de bons resultats que dans des situations 
favorables, comme des petites ozones ou le rayonnement 
diffuse est faible, ou des objets au contenu 
relativement homogene. II n' existe aucun procede simple 
qui permette par exemple de corriger le rayonnement 
diffuse a travers le thorax ou d'autres grandes zones 
anatomiques, dont l'examen est frequent mais qui sont 
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defavorables pour corriger le rayonnement diffuse en 
raison de leur volume meme et de 1 'heterogeneite due a 
la presence d'une structure d'os complexe et dont la 
capacite d' attenuation du rayonnement est tres 
differente de celle des tissus mous . 

Mentionnons enfin le brevet americain 
6 018 565 pour 1' expose d'une methode mixte, a « beam 
stop » et convolution. 

Un objet essentiel de 1' invention est de 
proposer un precede d' estimation et de correction de 
rayonnement diffuse qui puisse convenir pour des 
situations difficiles de radiographic . 

Le procede conforme a 1' invention est, sous 
sa forme la plus generale, un procede d'estimation d'un 
rayonnement diffuse provenant d'un rayonnement initial 
ayant traverse un objet en subissant une attenuation 
laissant passer un rayonnement total de mesure, 

caracterise par : 

- une prise d'une table de mesures d'un 
rayonnement diffuse, obtenue en faisant passer le 
rayonnement initial par un simulacre de 1' objet, 

un calcul de coefficients de 
transposition entre le simulacre et l'objet, d'apres le 
rayonnement initial, le rayonnement total de mesure a 
travers l'objet et un rayonnement total de mesure a 

travers le simulacre, v 

- et une ponderation de la table de mesures 
avec les coefficients de transposition. 

Avantageusement, le simulacre sera un bloc 
d'epaisseur constante et en une matiere homogene, ayant 
une attenuation semblable a une matiere de base de 
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l'objet ; en general la prise de table de mesure sera 
une selection dans une serie de tables de mesures de 
rayonnement diffuse, obtenues auparavant en faisant 
successivement passer le rayonnement initial a travers 
5 une serie respective de simulacres de l'objet, 
d'6paisseurs differentes mais constante ; et la 
selection sera faite par comparaison d'une valeur du 
rayonnement total de mesure a travers l'objet et d'une 
valeur du rayonnement total de mesure k travers les 
10 simulacres. 

Les coefficients de ponderation sont 
generalement des rapports de valeurs d'une meme 
f onctionnelle calculee pour l'objet et pour le 
simulacre. La f onctionnelle utilisee peut Stre 6gale au 

15 produit du rayonnement total de mesure par le 
logarithme du rapport de rayonnement total de mesure et 
du rayonnement initial. 

L' invention sera maintenant decrite en 
reference aux figures, parmi lesquelles : 

20 - la figure 1 est une vue generale d'une 

acquisition des mesures ; 

- la figure 2 est une vue d'une acquisition 
de calibration ; 

- et la figure 3 illustre les etapes du 

25 procede. 

Repor tons -nous d' aboard a la figure 1, oil un 
tube 1 de rayons X emet un faisceau 2 conique vers un 
objet 3 a examiner (ici un patient etendu sur une table 
4) puis, h travers lui, vers un reseau 5 plan de 
3 0 detecteurs 6 disposes en matrice. Les detecteurs 6 sont 
relics a un appareil d ' acquisition 7 et mesurent un 
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rayonnement diffuse qui se superpose au rayonnement 
primaire, seul convenable pour la radiographie • 

L' estimation du rayonnement diffuse a 
travers le patient 3 consiste tout d'abord a obtenir 
des tables bidimensionnelles ou nappes de rayonnement 
diffuse obtenues dans des circonstances comparables . 
Pour cela, on effectue des irradiations d'etalonnage a 
travers des simulacres 8 de l'objet 30 a examiner, 
conformement k la figure 2 : les conditions 
d' irradiation restent les m§mes, c'est-&-dire qu'on 
continue d'utiliser le tube 1, le faisceau 2, le reseau 
5 de detecteurs 6 et l'appareil d' acquisition 7, le 
simulacre 8 remplagant cependant le patient ; on a 
aussi ajoute une grille 9 de billes 10 de plomb entre 
le simulacre 8 et le reseau 5. II resulte de cette 
disposition que les rayons 11 passant par les billes 10 
sont compl£tement absorbes et que les regions 12 du 
reseau 5 situees dans le prolongement de ces rayons 11 
ont des detecteurs 6 qui ne mesurent que le rayonnement 
diffuse a ces endroits . II suffit de relever ,, y ces 
valeurs rnesurees et d'interpoler entre les regions 12 
pour estimer convenablement le rayonnement diffuse issu 
du simulacre 8 pour tous les detecteurs 6 du reseau 5 . 

Le simulacre 8 devrait §tre semblable a 
l'objet afin que les rayonnements diffuses par eux 
fussent identiques. Une similitude parfaite n'est pas 
realisable, et c'est pourquoi on se contente d'un 
simulacre 8 ressemblant a l'objet 3 et dont la nappe 
associee de rayonnement diffuse sera corrigee 
ulterieurement pour evaluer celle de l'objet. En 
pratique, le simulacre 8 peut etre un bloc d'une 
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matifere homogene et qui pr^sente le meme coefficient 
d' attenuation que la matiere de base de l'objet 3 : 
dans le cas d'un corps humain, compost pour l'essentiel 
de tissu mou, on sait que le plexiglas 
(polym^thacrylate) convient . 

Afin de permettre des mesures variees, on 
disposera en realite de plusieurs nappes de rayonnement 
diffuse, obtenues pour autant de simulacres 8, qui ne 
differeront que par leur epaisseur et done par la 
longueur du trajet parcouru par les rayons 11. Ces 
nappes seront enregistr^es dans une base de donnees 
pr<§alablement aux mesures utiles sur les objets 3 
radiographics, pour prendre une nappe de rayonnement 
diffuse comparable a celle d'un objet 3, on 
selectionnera en pratique une des nappes de la base de 
donn£es ou, mieux, une nappe qu'on aura obtenue par des 
calculs d' interpolation entre deux de ces nappes. Le 
critere de selection pourra §tre defini au moyen d'un 
rayon particulier 13 aboutissant a une region 14 du 
reseau 5 et qui ne passera ni par les absorbeurs 10 de 
la figure 2, ni par des tissus osseux du patient (ou 
plus g^neralement des portions de l'objet 3 dont les 
proprietes d' absorption sont differentes du materiau du 
simulacre 8) a la figure 1. Le rayonnement total, 
primaire et diffuse, regu par la region 14 apres avoir 
traverse chaque simulacre 8 servira d' index a la table 
de rayonnement diffuse correspondante, et la table 
selectionnee aura 1' index a une valeur identique au 
rayonnement total mesure a la region 14 a travers 
l'objet 3. Tout cela correspond au passage de 1 ' etat El 
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a l'etat E2 dans 1 ' or gani gramme de la figure 3, qu'on 
commence a commenter. 

La suite du procede consiste 

essentiellement en la correction de la table du 
5 rayonnement dif fusee ainsi selectionne pour 1'ajuster 
au mieux qu'on puisse esperer a la nappe de rayonnement 
reellement diffuse par 1'objet 3. Pour cela, on se sert 
de toutes les informations disponibles, c'est-a-dire du 
rayonnement total regu par les detecteurs 6 au-del^ de 
10 1'objet 3 comme du simulacre 8 selectionne. Ce 
rayonnement total etant note <I>t, le rayonnement diffuse 
<E>d, le rayonnement initial issu du tube 1 Oo et le 
rayonnement primaire 3>, la relation <t>t = O + Od est 
respectee . 

15 On est alors aux etats E3 et E4 de 

1 ' organigramme de la figure 3. Ensuite, on transforme 
les valeurs des rayonnements totaux Ot mesurees pour 
1'objet 3 et le simulacre 8 selectionne en leur 
appliquant des f onctionnelles . Plus precisement, il est 

20 connu dans 1 ' art que <I>d est proportionnel ci 
01og(0/<I>o) ; cette relation, qui est appelee la loi de 
Klein et Nishina, donne une allure generale du 
rayonnement diffuse, it defaut de son intensite. 

Le rayonnement initial O0 est connu ; le 

2 5 rayonnement primaire 4> ne 1'eSt pas, mais on consent a 
appliquer cette relation de fagon approchee en le 
remplagant par le rayonnement total <X>t^ c'est-a-dire 
que la f onctionnelle employee associe a chaque valeur 
mesuree du rayonnement total <I>t la valeur calculee 

30 Otlog(Ot/<E>o) , supposee proche du rayonnement diffuse 
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<I>d a cet endroit ; on est parvenu aux etats E5 et E6 de 
1 ' organigramme . 

L'etape suivante consiste a faire, pour 
chacun des detecteurs 6, le rapport des valeurs donnees 
par la f onctionnelle pour l'objet 3 et le simulacre 8 
selectionne selon i a formule 



K = — 1 °g(^ t / Oobbj et 



<E>t log(<I>t / Oo jsimulacre * Les coefficients de 
ponderation K ainsi obtenus serviront a deformer la 
nappe de rayonnement diffuse selectionnee a l'etat E2 
afin d'es timer celle de l'objet 3. Les resultats 
constituent encore une table bidimensionnelle ou une 
matrice ayant des dimensions identiques a celle des 
tables de rayonnement puisqu'elle est associee au 
reseau 5 de detecteurs 6. II est done possible et 
avantageux d'effectuer un filtrage numerique spatial de 
cette matrice en appliquant un filtre passe-bas qui 
corrige les coefficients K en ne conservant que les 
frequences les plus basses de leur variation et de les 
rendre ainsi probablement plus conformes a la realite 
puisque le rayonnement diffuse varie assez lentement 
d' un point a un autre. 

Quand la table des coefficients de 
ponderation definitifs, notes K' , a ete obtenue (a 
l'etat E7), elle sert a ponderer la table de 
rayonnement diffuse selectionnee auparavant a l'etat 
E2, pour obtenir une table de rayonnement diffuse par 
l'objet 3 (etat E8, qui constitue 1 ' estimation 
recherchee) ; la formule appliquee est 

<E>d objet=K' Od simulacre, Ces valeurs estimees Od 
objet pourront alors etre soustraites du rayonnement 
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total Ot mesur6 par les detecteurs 6 pour estimer le 
rayonnement primaire 3> et obtenir une image 
radiographique plus precise de l'objet 3. 

Ce procede s' applique aux radiographies a 
5 energie d' irradiation simple ou multiple ; dans le 
second cas, il est r6p6t6 separement pour chacune des 
energies employees . 

La f onctionnelle proposee ici n'est pas la 
seule qu'on puisse employer, et la f onctionnelle plus 
10 simple Od=k<I> (approch^e ici encore en Od=k<J>t) , k £tant 
une constante, pourrait aussi donner de bons r^sultats 
pour estimer Od. 

15 
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REVINDICATIONS 

1. Procede d' estimation d'un rayonnement 
diffuse provenant d'un rayonnement initial ayant 
traverse un objet (3) en subissant une attenuation 
laissant passer un rayonnement total de mesure, 
caracterise par : 

- une prise d'une table de mesures d'un 
rayonnement diffuse, obtenue en faisant passer le 
rayonnement initial par un simulacre (8) de 1' objet, 

un calcul de coefficients (K ' ) de 
transposition entre le simulacre et 1' objet, d'apres le 
rayonnement initial ($o) , le rayonnement total de 
mesure a travers 1' objet (<Dt objet) et un rayonnement 
total de mesure a travers le simulacre (Ot simulacre) , 

- et une ponderation de la table de mesures 
avec les coefficients de transposition. 

2. Procede d' estimation d'un rayonnement 
diffuse selon la revendication 1, caracterise en ce que 
le simulacre (8) est un bloc d'epaisseur constante et 
en une matiere homogene, ayant une attenuation 
semblable a une matiere de base de 1' objet. 

3. Procede d'estimation d'un rayonnement 
diffuse selon la revendication 1, caracterise en ce que 
la prise de table de mesures est une selection dans une 
serie de tables de mesures de rayonnement diffuse, 
obtenues en faisant success ivement passer le 
rayonnement initial a travers une serie respective de 
simulacres de 1' objet, qui sont des blocs d'epaisseurs 
differentes mais constante et en une matiere homogene, 
ayant une attenuation semblable a une matiere de base 
de 1' objet. 
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4. Procede d' estimation d'un rayonnement 
diffuse selon la revendication 3, caracterise en ce que 
la selection comprend une interpolation entre deux des 

tables de mesures . 

5. Procede d'estimation d'un rayonnement 
diffuse selon l'une quelconque des revendications 3 et 
4, caracterise en ce que la selection est faite par 
comparaison d'une valeur du rayonnement total de mesure 
a travers l'objet et d'une valeur du rayonnement total 
de mesure a travers les simulacres . 

6. Procede d'estimation d'un rayonnement 
diffuse selon la revendication 5, caracterise en ce que 
la comparaison est faite pour des rayons identiques 
(13) du rayonnement initial a travers l'objet et les 
simulacres, ne traversant que la matiere de base de 
1 ' ob j et . 

7. Procede d'estimation d'un rayonnement 
diffuse selon l'une quelconque des revendications 1 a 
6, caracterise en ce que les coefficients de 
ponderation sont des rapports de f onctionnelle 
identiques calculees pour l'objet et pour le simulacre. 

8. Procede d'estimation d'un rayonnement 
diffuse selon la revendication 7, caracterise en ce que 
les fonctionnelles sont egales au produit du 
rayonnement total de mesure par le logarithme du 
rapport du rayonnement total de mesure et du 
rayonnement initial . 

9. Procede d'estimation d'un rayonnement 
diffuse selon l'une quelconque des revendications 1 a 
8, caracterise en ce qu'il comprend une etape de 
filtrage passe-bas des coefficients de transposition, 
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arranges en une table superposable a la table de 
mesures . 

10. Procede de radiographic comprenant une 
etape de correction de mesures de radiographic par une 
sous tract ion d'un rayonnement diffuse estime selon le 
precede selon l'une quelcongue des revendications 
pr6c£dentes . 
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